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Abstract 

Recently, several studies of non-Hermitian Hamiltonians having VT symmetry have 
been conducted. Most striking about these complex Hamiltonians is how closely their 
properties resemble those of conventional Hermitian Hamiltonians. This paper presents 
further evidence of the similarity of these Hamiltonians to Hermitian Hamiltonians by 
examining the summation of the divergent weak-coupling perturbation series for the 
ground-state energy of the 'PT-symmetric Hamiltonian H = p 2 + jx 2 + i\x :i recently 
studied by Bender and Dunne. For this purpose the first 193 (nonzero) coefficients 
of the Rayleigh-Schrodinger perturbation series in powers of A 2 for the ground-state 
energy were calculated. Pade-summation and Pade-prediction techniques recently de- 
scribed by Weniger are applied to this perturbation series. The qualitative features of 
the results obtained in this way are indistinguishable from those obtained in the case 
of the perturbation series for the quartic anharmonic oscillator, which is known to be 
a Stieltjes series. 



1 Introduction 

Hamiltonians describing fundamental interactions traditionally possess two symmetries, the 
continuous symmetry of the proper Lorentz group and the discrete symmetry of Hermiticity. 
Lorentz invariance is a physical requirement. Hermiticity is a useful mathematical constraint 
that guarantees that the spectrum is real, although recent work shows that Hermiticity is 
only a sufficient condition and is not necessary for the reality of eigenvalues. From the 
assumptions of Lorentz invariance and positivity of the spectrum of the Hamiltonian one 
can prove the VCT theorem and thereby establish the physical symmetry of VCT invariance. 

What happens if we impose only the more physical symmetries of Lorentz invariance and 
VCT invariance when we construct a Hamiltonian? The constraint of VCT invariance is 
weaker than Hermiticity, so Hamiltonians having this property need not be Hermitian. While 
it has not been proved, there is compelling analytical and numerical evidence supporting the 
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conjecture that, except when PCT symmetry is spontaneously broken, the energy levels of 
such Hamiltonians are all real and positive The reality and positivity of the spectrum 

is apparently a consequence of the PCT symmetry of H. 

Many examples of PCT-symmetric Hamiltonians in quantum field theory have been 
studied In quantum mechanics, where the C operator is unity, many examples of 

PT-symmetric Hamiltonians have also been studied [p|-^9|. A simple example of such a 
quantum- mechanical Hamiltonian is H = p 2 + ix 3 . Hamiltonians such as this may be 
regarded as complex deformations of conventional Hermitian Hamiltonians. To understand 
this deformation we consider the Hamiltonian 

H = p 2 ~ (ix) 2+t , 

where e > 0. When e = 0, we have the conventional harmonic oscillator Hamiltonian, whose 
spectrum is real and positive. As e increases from 0, the entire spectrum of the Hamiltonian 
smoothly deforms as a function of e and remains real and positive for all positive values 
of e. Thus, these theories are in effect the analytic continuation of conventional quantum 
mechanics into the complex plane. 

These non-Hermitian theories exhibit some remarkable properties. Most interesting is 
that the expectation value of the operator x in quantum mechanics (and of the field (j> in 
quantum field theory) is nonzero when e > 0. This is true even for the p 2 — x 4 Hamiltonian 
that one obtains at e = 2 and it is also true for the —g<t> A scalar quantum field theory. The 
—g<j) A quantum field theory is particularly surprising because it has a positive real spectrum 
and exhibits a nonzero value of (<f>). In four-dimensional space-time it has a dimensionless 
coupling constant, is renormalizable, and is asymptotically free (and thus nontrivial). It 
may thus provide a useful setting to describe the Higgs particle 0. 

We are struck by the close similarity between the properties of non-Hermitian PT- 
symmetric quantum-mechanical Hamiltonians and conventional Hermitian Hamiltonians. 
Moreover, in mathematical terms, we are struck by the strong resemblance between self- 
adjoint Sturm-Liouville problems and these new complex Sturm-Liouville problems. The 
purpose of this paper is to present further evidence of this strong similarity by investigat- 
ing various aspects of Pade summation and Pade prediction of the Rayleigh-Schrodinger 
perturbation series for the ground-state energy of the complex 'PT-symmetric Hamiltonian 

H(X) = p 2 + ^x 2 + i\x 3 . (1) 

Note that this Hamiltonian is PT symmetric because under parity reflection V : p — > — p 
and P : x — > — x and under time reversal, which is an antiunitary operation, T : p — > — p, 
T : x — > x, and T : i — > —i. 

The large-order behavior of the divergent Rayleigh-Schrodinger perturbation series 

1 00 

EoW ~ 2 + 2 K)?n ( A ^ 0+ )< ( 2 ) 

n=l 

for the ground-state energy eigenvalue of the Hamiltonian in (Q) has already been examined 
in p(J, where the first 46 terms of the perturbation expansion had been generated using 
recursion formulas. It was observed there that the coefficients b n are all integers, that they 
alternate in sign, and that their magnitude grows rapidly with n. The first 10 coefficients 
are listed in Table [l]. 

To calculate the coefficients b n we make the ansatz that the wave function is a formal 
series in powers of the coupling constant A and that the coefficient of A" has the form of 
a Gaussian exp(— x 2 /A) times a polynomial of degree 3n in the variable x. The eigenvalue 
Eq(X) automatically appears as a series in powers of A 2 . Thus, for each additional coefficient 
in the series for E$ it is necessary to calculate two orders in powers of A for the wave function. 



On the Stieltjes nature of a "PT-symmetric perturbation expansion 



3 



Table 1: The first 10 coefficients b n in the perturbation expansion (g) for the 
ground state energy of the complex PT-symmetric Hamiltonian ([!]) . 



n 


b„ 


1 


11 


2 


-930 


3 


158836 


4 


-38501610 


5 


11777967516 


6 


-4300048271460 


7 


1815215203378344 


8 


-868277986898581530 


9 


464025598165231889260 


10 


-274145574452876905074540 



In [£0| it was pointed out that the Hamiltonian ([l]) describes a + 1 dimensional </> 3 field 
theory and that <j) 3 theories were the first quantum field theories in which the divergences 
of perturbation theory were studied [gjj. Using the standard methods for determining the 
large-order behavior of perturbation theory 22 , 23] it can be shown that the leading large- n 



behavior of the coefficients b n is given by the formula 

b n = (-l)™ +1 60 n+1/2 (27r)- 3/2 I> + l/2) [l + 0(l/n)] (n-»oo). (3) 



This asymptotic behavior was verified numerically in |20|. There, it was also shown that 
the first correction term to this leading asymptotic behavior, which is proportional to 1/n, 
is negative. Although divergent, the series in (||) is Borel summablc p4|-p7jl. If the factor 
of i were absent from the Hamiltonian (|l|), then the perturbation coefficients b n would not 
alternate in sign and the perturbation series would not be Borel summable. 

It is interesting that to derive the asymptotic formula in (|^) one must use dispersion- 
relation techniques that rely on crucial assumptions about the analyticity of the function 
Eq(X). These assumptions are justified for the Hermitian Hamiltonian of the anharmonic 
oscillator gH], 

Hifl) = P 2 + x 2 + /fe 4 . (4) 

However, the validity of these assumptions is unproved for the non- "PT-symmetric Ax 3 
oscillator. 

The eigenvalues £{0) of the quartic anharmonic oscillator possess several other properties 
which closely resemble those of the eigenvalues E(X) of the PT-symmetric Hamiltonian 
(Q). For example, the ground-state energy eigenvalue £o(P) of the quartic anharmonic 
oscillator possesses a divergent weak-coupling perturbation expansion, which also diverges 
factorially |2|,|§: 

OO 

£b(/9) ~ l + ^£„/3" 03-O+), (5) 

71=1 

B n = (-^^(z) r(n + l/2)[l + 0(l/n)] (n - oo) . (6) 
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A comparison of the large-order asymptotics (|§|) and (|^) shows that the two perturbation 
expansions @) and (||) possess the same rate of divergence if we choose 

/3 = 40 A 2 . (7) 



In view of these striking similarities between the ground-state eigenvalues Eq(X) and 
£o(P) it should be interesting to investigate what other similarities do exist. In particu- 
lar, we are interested in similarities that could provide evidence that the Pade summation 
of the divergent perturbation series (|^) for the ground-state energy of the "PT-symmetric 
Hamiltonian ([]]) converges. 

In the case of the ground-state energy shift A£o((3) of the quartic anharmonic oscillator, 
which is defined by 



£ o (0) = 1 + 0A£b(/3), (8) 

it was shown rigorously by Simon (Theorem IV. 2.1 of |^8| ) that the corresponding pertur- 
bation series is a Stieltjes series. This has some far-reaching consequences. In the case 
of Stieltjes series, Pade approximants possess a highly developed convergence theory, as we 
discuss in detail in Section In particular, the Stieltjes nature of the perturbation series (^) 
guarantees that certain subsequences of the Pade table converge to a uniquely determined 
Stieltjes function. 

Although we cannot prove it rigorously, we believe that for the 'PT-symmetric Hamilto- 
nian (|]) the ground-state energy shift A£^ (A 2 ) defined by 

So (A) = i + A 2 A£ (A 2 ), (9) 

in which E is considered as a function of A 2 , is also a Stieltjes function. This implies that 
the corresponding perturbation series is a Stieltjes series. 

It is the intention of this article to provide numerical evidence supporting this conjecture. 
We do this by comparing Pade summations and Pade predictions (see [ |3l| and references 
therein) of the perturbation expansions for the ground-state energy shift Ai?o(A 2 ) and for 
the analogous ground-state energy shift A£q([3) of the quartic anharmonic oscillator. 

Here, one might argue that one should also investigate the summation of the pertur- 
bation expansion for the ground-state energy shift Ai?o(A 2 ) with the help of the sequence 
transformations that were described in Section 7 and 8 of |32] and which produced very 
good results in the case of the anharmonic oscillators |50|, [33], £34| . However, the convergence 
theory of these sequence transformations, which in the case of power series also produce 
rational approximants, is still very much in its infancy and no theoretical results concerning 
the transformation of Stieltjes series are known so far. Consequently, we would only produce 
numbers without gaining any further mathematical insight. 

In Section |[ we present the relevant details of Pade approximants, the computation of 
Pade approximants by means of Wynn's recursive epsilon algorithm |35f| , and we discuss 
Stieltjes series and their associated Stieltjes functions. In Section [|, we show that the Pade 
summation of the perturbation expansions for Ai?o(A 2 ) and A£q{0), respectively, produces 
results of identical quality if the two coupling constants A and /3 satisfy (0). In Section 
||, we discuss the prediction of unknown perturbation coefficients with the help of Wynn's 
epsilon algorithm, and we show that the coefficients of the perturbation expansions for 
either AE (X 2 ) or A£ (/3) can be predicted equally well. Finally, in Section || we give a 
brief summary. 
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2 Pade Approximants, Wynn's Epsilon Algorithm, and 
Stieltjes Series 

In recent years, Pade approximants have become the standard tool in theoretical physics to 
overcome problems with slowly convergent or divergent power series. Accordingly, there is 
a vast literature on the mathematical properties of Pade approximants as well as on their 
applications in theoretical physics. Any attempt to provide a reasonably complete biblio- 
graphy would be beyond the scope of this article (see for example the extensive bibliography 
compiled by Brezinski |3(|). We just mention that the popularity of Pade approximants in 
theoretical physics can be traced back to a review by Baker j37j, that the first specialized 
monograph on Pade approximants is due to Baker |3S| ] , and that currently the most complete 
source of information is the second edition of the monograph by Baker and Graves-Morris 
J39| . In addition to treatments in more mathematically oriented books on continued fractions 
and related topics jh)] Q, Pade approximants are also discussed in book on mathematical 
and theoretical physics, for example in Section 8 of the book by Bender and Orszag |2(i|] , or 
in Part III of a book by Baker on critical phenomena p^ j. Then, there is a book by Pozzi \il ] 
on the use of Pade approximants in fluid dynamics. Finally, there is even a monograph [46] 
and two articles |47],fl8) on the history of Pade approximants and related topics. 

A Pade approximant P} n (z) to a function / possessing a (formal) power series expansion 

oc 

/(*) = £7**", (io) 

!/=0 

which may converge or diverge, is the ratio of two polynomials Ai(z) and B m (z) of degrees 
I and to in z (p. 383 of ^(|): 

j . . Ai(z) a + a\z + a 2 z 2 + ... aiz 1 

m[Z) ~ B m (z) ~ l + b 1 z + b 2 z2 + ...b m z™- [ ' 

An alternative notation for Pade approximants, which is used in the books by Baker and 
Graves-Morris |38|,|3£]], is Pt n (z) — [l/m]f(z). This notation is usually simplified further to 
Pm(z) — [//to] if explicit references to / and z are not necessary. 

The coefficients a\ and of the polynomials Ai(z) and B m (z) in jll| ) are chosen in 
such a way that the Taylor expansion of f(z) and of its Pade approximant agree as far as 
possible: 

f(z)-Pi(z) = 0(z i+m+1 ) (*-0). (12) 

This asymptotic error estimate leads to a system of linear equations by means of which the 
coefficients ao, a%, ... , ai and b%, 62, • • • , b m in ( |i~T| ) can be computed [p8|p^| . Moreover, 
several algorithms are known which permit a recursive computation of Pade approximants. 
A discussion of the merits and weaknesses of the various computational schemes can for 



instance be found in Section II. 3 of the book by Cuyt and Wuytack |49|. 

Probably, the best known recursive algorithm for Pade approximants is Wynn's epsilon 
algorithm pq |: 

£ { "l = , 4 n) = s « . « e N , (13a) 
4" + } i = 4-i 1) + l/[4 n+1) -4" ) ]. k,ne N . (13b) 

A compact FORTRAN program for the epsilon algorithm as well as the underlying compu- 
tational algorithm is described in Section 4.3 of |32]| . 
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If the input data £q"^ = s n of Wynn's epsilon algorithm are the partial sums 



n 

fn(z) = (14) 

of the formal power series ( |Io|) according to = f n (z), then the elements e^b with even 
subscripts are Pade approximants to f(z) according to |35| 

4? = Pk +n (z)- (15) 

(n) 

The elements e^fe+i with odd subscripts are only auxiliary quantities, which diverge if the 
whole transformation process converges. 

The epsilon algorithm is a useful numerical algorithm that is applied successfully in a 
large variety of different fields. Accordingly, there is an extensive literature dealing with it. 
A fairly complete coverage of the older literature can be found in a book by Brezinski |30| . 
It may be interesting to note that the epsilon algorithm is not restricted to scalar sequences 
but can be generalized to cover vector sequences. A recent review of these developments 
can be found in jsj]] . 

If one tries to sum a divergent power series by converting its partial sums (|l4| ) to Pade 
approximants, it is usually a good idea to use diagonal Pade approximants, whose numerator 
and denominator polynomials have equal degrees. If this is not possible one should use Pade 
approximants with numerator and denominator polynomials whose degrees differ as little 
as possible. If we use the epsilon algorithm for the computation of the Pade approximants, 
then (|l5| ) implies that we should use the elements of following staircase sequence in the Pade 
table as approximations to f(z) (see Eq. (4.3-7) of |^| ): 

Pg(z),Pfa),P}(z),... ,P»{z),P» +1 {z),P»+t{z),... . (16) 

This staircase sequence exploits the available information optimally if the partial sums f n (z) 
with n > are computed successively and if, after the computation of each new partial sum, 
the element of the epsilon table with the highest possible even subscript is computed. With 
the help of the notation [xj for the integral part of x, this staircase sequence can be written 
compactly as follows: 

(n-2[n/2J) _ p n-|[n/2]|/A fn — d 1 9 ^ (-[7\ 

£ 2[[n/2l — r ln/2\ \ Z ) {Tl — V, 1, /,,...) . [II) 

As remarked above, Pade approximants are now used almost routinely to overcome 
problems with slowly convergent or divergent power series. Bence, their practical usefulness 
is beyond question. Bowever, from a theoretical point of view, the situation is not so 
good. So far, a completely satisfactory general convergence theory of Pade approximants 
for essentially arbitrary power series does not exist. 

Nevertheless, there is a special class of series, the so-called Stieltjes series, which possess a 
highly developed and elegant convergence theory. In this section we will only discuss those 
properties of Stieltjes series and Stieltjes functions that are needed to provide numerical 
evidence that the perturbation expansion for the ground-state energy shift AEq(X 2 ) of the 
complex 'PT-symmetric Bamiltonian considered as a function of A 2 , is a Stieltjes series. 
Detailed discussions of the properties of Stieltjes series and their special role in the theory 
of Pade approximants can be found in Section 8.6 of |^6| or in Section 5 of ]39[ ]. 

A function F{z) with z G C is called a Stieltjes function if it can be expressed as a 
Stieltjes integral according to 

r°° d$(t) 

F(z) = j o ^ [|arg(*)|<7r]. (18) 
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Here, $(£) is a bounded, nondecreasing function taking infinitely many different values on 
< t < oo. Moreover, the moment integrals 



f l d$(i) (neN ) 



(19) 



must be real and finite for all finite values of n. A Stieltjes function can be expressed by its 
corresponding Stieltjes series: 



F{z) = £(-1)"/^ 



(20) 



Whether this series converges or diverges depends on the behavior of the Stieltjes moments 
li n as n — > oo. 

In a typical Stieltjes summation problem, as it occurs in the context of divergent per- 
turbation expansion, only the numerical values of a finite number of Stieltjes moments /i„ 
are known. Thus, one has to find a way of constructing an approximation to the unknown 
Stieltjes function F{z) from a finite string of moments. 

Of course, one would also like to have some theoretical evidence that F(z) exists and 
is uniquely determined by the Stieltjes moments {/in}^Lo- Many necessary and sufficient 
conditions that guarantee this are known in literature. Unfortunately, it is by no means easy 
to apply them. For example, a necessary condition that the series ( |20| ) is indeed a Stieltjes 
series is that the Hankel determinants 



D(m, n) 



■1 Mm+2 



Mm+n+1 
Atm+2n 



(21) 



are positive for all m,n > (see Theorem 5.1.2 on p. 197 of p9[). 

Then, there is a sufficient criterion, the so-called Carleman condition (see p. 410 of j2(| 
or pp. 239 - 240 of p9|]), which requires that the series X^i(Mj) -1 diverges and thus 
limits the admissible growth of the moments \i n as n — > oo ]52] ]. If the condition D(m, n) > 
on the determinants as well as the Carleman condition are both satisfied, then the Pade 
approximants P™ +:, (z) constructed from the partial sums of the moment expansion ( po| ) 
converge for every j > —I to the corresponding Stieltjes function F(z) as m — > oo (see for 
example Theorem 5.5.1 on p. 240 of |p9[ ). 

If only a finite number of moments are known, it is impossible to prove that D(m, n) > 
holds for all m,n > 0, and it is also not possible to prove rigorously that the Carleman 
condition is satisfied, although we would like to emphasize that the large-order formula (^), 
which was verified numerically in p0| , is in agreement with the Carleman condition. 

Therefore, we prefer an indirect approach in order to provide evidence that the pertur- 
bation series for the energy shift AEq(\ 2 ) of the PT-symmetric Hamiltonian (|l) is indeed 
a Stieltjes series. For that purpose, let us assume that the moment expansion (20(), whose 
Stieltjes nature we want to establish, is a Stieltjes series. Pade approximants to Stieltjes 
series possess a highly developed convergence theory, and many conditions and inequalities 
are known that Pade approximants to a Stieltjes series must satisfy. For example, Pade 
approximants constructed from the partial sums 



i/=0 



(22) 
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of the moment expansion (|20j) for a Stieltjes function F(z) satisfy for z > the following 
inequalities (Theorem 15.2 on p. 215 of |3q|): 

{P™# +1 (z) - P™ +j (z)} > 0, (23) 

(-iy +1 [P% +j {z) - P%±?\z)} > 0, (24) 

P%(z) > F(z) > P^iz). (25) 

It follows from inequality ( p3|) that the Pade sequence {P™ + i (z)}°°_ Q is increasing for 
z > if j is odd, and it is decreasing if j is even. Moreover, if we set j = — 1 in (|^) and 
replace m by m + 1, we obtain the inequality 

P™ +1 (z) > P£ +1 (z) (me No). (26) 

Thus, if we use Wynn's epsilon algorithm ( jllf ) to convert the partial sums (g2) to Pade 
approximants and choose the approximants to F(z) according to jl^), then it follows from 
inequalities ( p5|) and ( |26|) t hat these Pade approximants satisfy the following inequality if 
the moment expansion (EG) for F(z) is a Stieltjes series: 

P™ +1 (z) < F(z) < P™+l(z) (z > 0, m e No) . (27) 

Thus, the approximants (JT^) produced by Wynn's epsilon algorithm yield for z > two 

nesting sequences P™ +1 (z) = e^m an d Pm+i( z ) = £ 2m+2 °f lo wer an d upper bounds to 
F(z) if the the moment expansion ( |20| ) is a Stieltjes series. 

If only the numerical values of a finite number of Stieltjes moments are available, then it 
is of course not possible to prove rigorously that the series under consideration is a Stieltjes 
series. Nevertheless, we can provide considerable evidence that this hypothesis is true if 
inequality (|27j ) is valid in all cases that can be checked. 



3 Summation Results 

In this section we want to show that the Pade summation of the perturbation expansion 

oo 

A£ (A 2 ) = J2 bu+i\ 2u (28) 

v=0 

for the ground-state energy shift of the of the T'T-symmetric Hamiltonian (0) and of the 
perturbation expansion 

oo 

ASo(p) = Y, Bv+iF ( 29 ) 

for the ground-state energy shift of the quartic anharmonic oscillator yield results of virtually 
identical quality if the two coupling constants A and (3 satisfy (0). Moreover, we want to 
demonstrate numerically that inequality (p7|), which is satisfied in the case of the Pade 
summation a Stieltjes series, is apparently also satisfied. In our summation calculations, we 
used all coefficients b u and £>„ with < v < 193 which had been computed recursively. 

In this article, we compute all Pade approximants with the help of Wynn's epsilon 
algorithm (||). Thus, the partial sums 

n 

!/=0 
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and 

n 

a n ((3) = B »+^ U ( 31 ) 

are used as input data for Wynn's epsilon algorithm according to = s n (\ 2 ) or = 
a n {S), respectively, and the approximations to the energy shifts are chosen according to 



Table 2: Pade summation of the perturbation expansions ( pq ) and ( p9[ ) for the ground-state 
energy shifts Ai? (A 2 ) and A£ (/3) with A = 1/7 and (3 = 40/49, respectively. 



n 


Sn(A 2 ) 




IW2] V 1 


[n/2] 


n 
u 


0.110 


10 002 


0.750 


10 ooo 


1 1 nn nnn nnn nnn nnn 
11. uu uuu uuu uuu uuu 


n 7^n nnn nnn nnn nnn 
u. / ou uuu uuu uuu uuu 


1 


-0.798 


10 ooi 


-0.321 


10 ooo 


-7.97 959 183 673 469 


-0.321 428 571 428 571 


2 


0.582 


10 002 


0.315 


10 ooi 


6.76 871 520 405 468 


0.497 075 017 205 781 


3 


-0.269 


10 003 


-0.133 


10 002 


3.14 452 476 154 168 


0.283 471 705 042 096 


4 


0.177 


10 uu 


0.861 


10 ' 


5.92 770 890 838 469 


0.444 962 648 249 413 


5 


-0.134 


10 005 


-0.639 


10 003 


4.84 920 642 167 536 


0.379 736 282 027 717 


50 


0.153 


10 072 


0.684 


10 070 


5.52 416 958 165 793 


0.419 249 574 461 710 


51 


-0.964 


10 073 


-0.432 


10 072 


5.52 416 451 428 038 


0.419 249 241 261 250 


52 


0.620 


10 075 


0.278 


10 074 


5.52 416 888 260 688 


0.419 249 527 748 761 


53 


-0.407 


10 077 


-0.182 


10 076 


5.52 416 531 636 255 


0.419 249 293 076 390 


54 


0.272 


10 079 


0.122 


10 078 


5.52 416 839 738 891 


0.419 249 495 310 895 


101 


-0.210 


1Q 172 


-0.939 


10 170 


5.52 416 721 141 847 


0.419 249 415 925 473 


102 


0.263 


10 m 


0.118 


1Q 173 


5.52 416 721 422 990 


0.419 249 416 112 202 


103 


-0.334 


10 176 


-0.149 


1Q 175 


5.52 416 721 178 460 


0.419 249 415 949 529 


104 


0.427 


10 178 


0.191 


10 177 


5.52 416 721 397 212 


0.419 249 416 094 862 


105 


-0.552 


10 180 


-0.247 


10 179 


5.52 416 721 206 667 


0.419 249 415 968 069 


150 


0.318 


10 279 


0.142 


10 278 


5.52 416 721 306 531 


0.419 249 416 033 824 


151 


-0.590 


10 281 


-0.264 


10 280 


5.52 416 721 305 477 


0.419 249 416 033 119 


152 


0.110 


10 284 


0.493 


10 282 


5.52 416 721 306 436 


0.419 249 416 033 760 


153 


-0.207 


10 286 


-0.928 


10 284 


5.52 416 721 305 579 


0.419 249 416 033 187 


154 


0.392 


10 288 


0.175 


10 287 


5.52 416 721 306 359 


0.419 249 416 033 708 


188 


0.371 


10 367 


0.166 


1Q 366 


5.52 416 721 306 035 


0.419 249 416 033 489 


189 


-0.860 


10 369 


-0.385 


10 368 


5.52 416 721 306 009 


0.419 249 416 033 472 


190 


0.201 


10 372 


0.898 


10 370 


5.52 416 721 306 033 


0.419 249 416 033 488 


191 


-0.471 


10 374 


-0.211 


10 373 


5.52 416 721 306 011 


0.419 249 416 033 474 


192 


0.111 


10 377 


0.496 


10 375 


5.52 416 721 306 031 


0.419 249 416 033 487 



In Table p| we present illustrative results of the Pade summation of the perturbation 
expansion (Eq) for the ground-state energy shift Ai?o(A 2 ) with A = 1/7 and of the pertur- 
bation expansion ( ps| ) for the ground-state energy shift A£q(P) with (3 = 40/49. Thus, the 
two coupling constants A and j3 satisfy (R). This implies that the two perturbation expan- 
sions, whose partial sums are displayed in columns 2 and 3, should show the same rate of 
divergence. 
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Here, we must remember that the larger-order estimates (|3|) and (|6|), respectively, imply 
that the partial sums ([30|) for the ground-state energy shift AEo(X 2 ) of the PT-symmetric 
Hamiltonian are for all n > at least one order of magnitude greater than the partial 
sums ([n]) for the ground-state energy shift A£q(/3) of the quartic anharmonic oscillator. 
Otherwise, the observed rates of divergence in columns 2 and 3 are virtually identical. 

Moreover, the Pade approximants in columns 4 and 5 apparently satisfy the inequality 
(p7|), which holds if the series to be transformed is a Stieltjes series. If the index n in column 
1 is even (n = 2m), then the diagonal Pade approximants P™(A 2 ) and P™(/3) provide upper 
bounds that strictly decrease with increasing m and, if n is odd (n = 2m + 1), then the 
Pade approximants P™ +1 (A 2 ) and P™ +1 {/3) provide lower bounds that strictly increase with 
increasing m. 

We have done analogous summation calculations also for many other values of the cou- 
pling constants A and (3. Of course, the performance of the Pade summations depend very 
much on the size of the coupling constants. For smaller values of A and the corresponding 
f3, convergence is better than in Table ||, whereas for A = 1 and f3 = 40 only the first digit 
of the summation results stabilize. For larger values of A and f3, Pade summation produces 
only relatively crude upper and lower bounds. However, to emphasize that the typical qual- 
itative features of the summation results in Table || - the same rate of divergence of two 
perturbation series and the occurrence of strictly decreasing upper bounds P™ and strictly 
increasing lower bounds P™ +1 ~~ "were consistently observed in all cases considered. 

Thus, Wynn's epsilon algorithm is apparently unable to detect any substantial difference 
between the perturbation series (^9| ) for A£o(/3), whose Stieltjes nature was established 
rigorously by Simon (Theorem IV. 2.1 of [§8|), and the perturbation series (^8|) for APo(A 2 ), 
whose Stieltjes nature we conjecture. 



4 Pade Predictions 

As shown by countless articles from all branches of physics, Pade approximants have become 
the standard tool to overcome problems with slowly convergent or divergent power series. 
However, Pade approximants have other useful features that are not as well known yet. 
For example, Pade approximants can be used to make predictions for higher-order series 
coefficients that were not used for the construction of the approximant. 

On a heuristic level the prediction capability of Pade approximants, which was apparently 
first noted and used by Gilewicz ]53| ], can be explained quite easily. Let us assume that 
the partial sums (|l4|) of the power series for some function f(z) are to be converted to 
Pade approximants. Then, the accuracy-through- order relationship ( |l2| ) implies that a Pade 
approximant P^z) to f(z) can be expressed as the partial sum fi +m (z) from which it was 
constructed plus a term z l+m+1 lZ l m (z) , which was generated by the transformation of the 
partial sum to the rational approximant: 

PL(?) = E + z l+m+1 K l m (z) = fl+m(z) + z l + m+1 K\ n {z). (32) 

i/=0 

Similarly, the power series (|io| ) can be expressed as follows: 

l-\-m oo 

f(z) = J2^ Z " + zl+m+1 E H+m+u+xZ v = fi +m (z) + z l+m+1 T l+m+1 {z) . (33) 

i/=0 w=Q 

Let us now assume that the indices / and m so large that the Pade approximant Pl n (z) 
provides a sufficiently accurate approximation to f(z). Then, the Pade transformation 
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term lZ l m (z) must also provide a sufficiently accurate approximation to the truncation error 
J-i+ m +i{z) of the power series. In general, we have no reason to assume that 

K l m {z) = Fl +m+ i{z) (34) 

might hold exactly for finite values of I and to. Consequently, Taylor expansions of TZ l m (z) 
and J-;+ m +i(z), respectively, will in general produce different results. Nevertheless, the 
leading coefficients of the Taylor expansion for TZ m {z) should in such a case provide suf- 
ficiently accurate approximations to the corresponding coefficients of the Taylor series for 

It is important to note that this prediction capability of Pade approximants does not 
depend on the convergence of the power series expansions for H l m {z) and Ti +m +i{z), respec- 
tively, which was used implicitly in our heuristic reasoning given above. Pade approximants 
are able to make predictions about series coefficients even if the power series (|l0|) for / as 
well as the power series expansions for 72.J™ and Fi+ m +i{z) are only asymptotic as z — > 0. 
This fact explains why the prediction capability of Pade approximants can be so very useful 
in the case of violently divergent perturbation expansions (see pl|, p4|-|72|] and references 
therein). 

Theoretically, very little is known about the prediction of series coefficients that were not 
used for the construction of the Pade approximant. A notable exception are again Stieltjes 
series for which inequalities are known. 

Let us assume that the partial sums ( |22] ) of the moment expansion for some Stieltjes 
function F(z) are to be converted to Pade approximants. Such a Pade approximant P^z) 
possesses the following power series expansion: 

oo 

pLW = £(-i)>l' /m1 ^. (35) 

!/=0 

In Theorem 5.2.7 on p. 220 of J39| it was shown that for all n G No and for all I > m — 1 
the coefficients in © are bounded in magnitude by the Stieltjes moments /x„ in ( |2C| ) 

according to 

< /4! /ml < M« • (36) 

This inequality can be used to analyze the Stieltjes nature of a moment expansion of 
the type of (^0|). With the help of computer algebra systems like Maple or Mathematica 
it is possible to construct Pade approximants Pm( z ) m an unspecified symbolic variable 
z, and this can even be done free of rounding errors if the coefficients of the series to be 
transformed are exact integers like the coefficients b n in (Q) or exact rational numbers like 
the coefficients B n in (||). In the next step, a the leading part of a power series expansion 
of the Pade approximant must be constructed. The resulting series coefficients can then be 
compared with the corresponding coefficients of the moment expansion. 

Again, this poses no principal problems for computer algebra systems like Maple and 
Mathematica. However, the accuracy-through-order relationship ([l2]) implies that inequality 
(|36|), which is to be checked, is by default satisfied for all indices n < I + to, and only for 
n > I + to + 1 we obtain useful information about the Stieltjes nature of the moment 
expansion, from which the Pade approximant was constructed. Thus, if I + m becomes 
large, the brute-force approach based on computer algebra systems becomes very demanding 
both with respect to computer time and memory because it requires both the symbolic 
construction of complicated Pade approximants and also symbolic differentiations of very 
high orders. 

These computational problems can be simplified considerably with the help of a recently 
derived recursive scheme (Section 3 of pl|) that permits a direct calculation of the transfor- 
mation term 7Z l m (z) in ( |32| ) if the corresponding Pade approximant P [ m {z) can be computed 
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with the help of Wynn's epsilon algorithm according to (pL 5|) ; that is, for Pade approximants 
of the type P k+n {z) with k, n G N . 

It follows from the accuracy-through-order relationship ( |l2| ) in combination with ( |l5| ) 
that e 2 ^ can be expressed as follows if the partial sums Q of the power series for some 
function f(z) are used as input data for Wynn's epsilon algorithm (fiT 



£ 2k - Jn+2k{Z) + Z lfi 2k (Z) 



(37) 



The quantities ip^ (z) can be computed with the help of the recursive scheme in Eq. (3.15) 
of pl| , which uses the coefficients j n of the power series (|l0|) for f(z) as input data, as 
follows: 



fo ( z 
(")/ 



(«) / 



(n) 



^2fe+2^ 



where 



= 0, neNo, 
[7^+2] 2 

Jn+1 - ln+2Z ' 

(«) 

(n+2)/ \ a 2fc- 

( z ) + 



2fe+2 



7n+2fe+l + fy^OO 



7n+2fe+2 + ^^™ +1) (z) 



fcefj, n G N , 



1 X («+!)/ A 
7» + 2fc+2+<V 2fc '(Z) 

7n+2fc+l + Z¥>2fc \ Z ) ~ ^2k-2 ( Z ) 



7n+2fe+l + 5(p$(z) 



(n+1)/ \ yn-r^j 1 \ 

7n+2fe+l + Zip\ k (Z) - tp y 2k _ 2 ' (Z) 



(n+2). 



(38a) 
(38b) 

(38c) 
(38d) 

(38e) 
(39) 



The rational function ip^i \z) can be expressed as a power series in z according to 



^2k\ Z ) 



<n,2k) , (n,2k) 

9o + 9i z 



(n,2k) 2 , 
92 Z + 



g (n,2k) z , 



(n,2k) 



(40) 



of this series expansion can be used to predict the coefficients 



The coefficients g, 

ln+2k+v+i with v > of the power series expansion ([L0|) for f[z) that were not used 
for the construction of either or ip^(z). 

(n) 

Thus, we can compute the rational function ip 2k with the help of the recursive scheme 
( |38| ) in the case of the perturbation series ( |28| ) for AEq(X 2 ) and ( p9| ) for A£q(/3). In the 
case of the "PT-symmetric perturbation series we obtain the expansion 

f2k\^) = b ( t 2k) + bi n,2fe) A 2 + b { ^ 2k) A 4 + . . . + A 2 " + . . . , (41) 

and in the case of the perturbation series for the quartic anharmonic oscillator we obtain 

. . + B { J l/2k) 13" + ... . (42) 



<4;)(/3) = B^ 2k) + B[ n ' 2k) [3 + 4"' 2fe) /3 2 



The coefficients bu and g(, n ' 2fc ) with v > can be used to obtain predictions for the 
coefficients 6 n +2fc+i/+2 and B n +2k+v+2, respectively, that were not used for the construction 
of the Pade approximants e 



{n) 

2k 



P k +n or the transformation terms ip^ 2k . 
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In Table |^ we compute the rational function ip^ 6 ' (A 2 ) corresponding to the diagonal Pade 
approximant = P33 (A 2 ) with the help of the recursive scheme ( |38| ) from the coefficients 
b v with 1 < v < 67. The resulting expansion coefficients fry 0,66 -* defined in (|d]) with v > 
provide predictions to the coefficients 6^+68 of the perturbation series ( |28| ) . All calculations 
for Table || were done free of rounding errors using the exact rational arithmetics of Maple. 
Only in the final step the coefficients were converted to floating point numbers for the sake 
of readability. 

Table 3: Predictions bn' 66 ^ for the coefficients b n+es of the perturbation series ( p8|) for 
AE (X 2 ) with < n < 15. 



n n On+68 






-0 


118 


625 


502 


281 


564 


111 


353 


10 217 


-0 


118 


625 


502 


281 


564 


111 


358 


10 217 


1 





487 


707 


952 


691 


623 


584 


397 


10 220 





487 


707 


952 


691 


623 


585 


158 


10 220 


2 


-0 


203 


437 


822 


070 


101 


216 


978 


10 224 


-0 


203 


437 


822 


070 


101 


222 


504 


10 224 


3 





860 


803 


267 


021 


875 


481 


138 


10 227 





860 


803 


267 


021 


875 


756 


369 


10 227 


4 


-0 


369 


393 


498 


548 


727 


222 


559 


10 231 


-0 


369 


393 


498 


548 


728 


279 


960 


10 231 


5 





160 


732 


212 


082 


002 


560 


522 


10 235 





160 


732 


212 


082 


005 


901 


209 


10 235 


G 


-0 


709 


026 


471 


212 


486 


114 


145 


10 238 


-0 


709 


026 


471 


212 


576 


489 


701 


10 238 


7 





317 


020 


667 


799 


578 


470 


271 


10 242 





317 


020 


667 


799 


793 


728 


631 


10 242 


8 


-0 


143 


648 


198 


373 


426 


854 


924 


10 246 


-0 


143 


648 


198 


373 


887 


496 


043 


1Q 246 


9 





659 


514 


281 


085 


804 


565 


498 


10 249 





659 


514 


281 


094 


798 


452 


336 


10 249 


10 


-0 


306 


750 


687 


264 


795 


900 


309 


1Q 253 


-0 


306 


750 


687 


281 


012 


588 


650 


10 253 


11 





144 


514 


693 


689 


642 


646 


364 


10 257 





144 


514 


693 


716 


909 


093 


737 


10 257 


12 


-0 


689 


498 


329 


409 


437 


387 


151 


10 260 


-0 


689 


498 


329 


840 


371 


816 


155 


10 260 


13 





333 


104 


548 


293 


054 


144 


923 


10 264 





333 


104 


548 


937 


521 


023 


558 


10 264 


14 


-0 


162 


924 


769 


352 


053 


020 


131 


10 268 


-0 


162 


924 


770 


269 


205 


895 


837 


10 268 


15 





806 


654 


532 


549 


091 


198 


441 


10 271 





806 


654 


545 


029 


445 


531 


410 


10 271 



The results in Table show that the expansion coefficients b v ' with v > provide 
already remarkably accurate predictions to the corresponding coefficients 6„+68- Moreover, 
the Stieltjes inequality (^) is satisfied in all cases. 

In Table ^ we do the same calculations as in Table but this time for the ground- 
state energy shift A£o(/3) of the quartic anharmonic oscillator. Thus, the expansion of the 
rational function ip 1 ^ (/3), which is computed from the coefficients B v with 1 < v < 67 of 

the perturbation series (|29|), provides predictions Bv'^ to the coefficients #„ + 68- 

A comparison of Tables || and ^ shows that their qualitative features are identical. In 
this context it is quite remarkable that although the coefficients b„ grow significantly more 
rapidly in magnitude than the coefficients £>„, which follows from the large-order estimates 
(||) and (^), Pade prediction nevertheless yields results of virtually identical quality. 

The Pade prediction of unknown series coefficients based on the recursive scheme ( |3^ ) 
is certainly computationally simpler than the straightforward symbolic computation and 
expansion of Pade approximants. The rational function ip^} (z) has a simpler structure than 
£ 2k ~ ^ > fc + "( z )i an d the first n + 2k symbolic differentiations can be avoided. Nevertheless, 
the recursive symbolic computation of the rational function (z) from the coefficients 
70, 71, ... , 7n+2fc of the power series ( |lo| ) can become quite demanding, in particular if k 
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Table 4: Predictions si°' 66) 
AE (/3) with < n < 15. 



for the coefficients £> n +66 of the perturbation series 



for 



B ( 0> 66) 



B 



M+6S 






-0 


243 


941 


384 


991 


118 


295 


771 


10 108 


-0 


243 


941 


384 


991 


118 


295 


782 


10 108 


1 





250 


725 


042 


695 


070 


353 


544 


10 no 





250 


725 


042 


695 


070 


353 


955 


10 no 


2 


-0 


261 


457 


030 


278 


874 


510 


535 


10 112 


-0 


261 


457 


030 


278 


874 


517 


978 


10 112 


3 





276 


569 


040 


522 


183 


341 


803 


10 114 





276 


569 


040 


522 


183 


434 


367 


10 U4 


4 


-0 


296 


701 


814 


375 


736 


021 


569 


10 116 


-0 


296 


701 


814 


375 


736 


909 


442 


10 U6 


5 





322 


749 


390 


515 


363 


534 


568 


10 118 





322 


749 


390 


515 


370 


538 


244 


10 U8 


G 


-0 


355 


923 


577 


678 


312 


610 


650 


10 120 


-0 


355 


923 


577 


678 


359 


918 


630 


10 12O 


7 





397 


845 


013 


388 


761 


856 


087 


10 122 





397 


845 


013 


389 


043 


208 


304 


10 122 


8 


-0 


450 


670 


140 


529 


734 


425 


361 


10 124 


-0 


450 


670 


140 


531 


237 


820 


183 


10 124 


9 





517 


267 


603 


130 


982 


724 


472 


10 126 





517 


267 


603 


138 


312 


496 


560 


10 126 


10 


-0 


601 


463 


530 


952 


366 


420 


452 


10 128 


-0 


601 


463 


530 


985 


369 


171 


801 


10 128 


11 





708 


383 


831 


448 


709 


420 


419 


10 130 





708 


383 


831 


587 


280 


706 


922 


10 130 


12 


-0 


844 


934 


259 


800 


460 


726 


512 


10 132 


-0 


844 


934 


260 


347 


380 


950 


714 


10 132 


13 





102 


047 


769 


191 


033 


928 


036 


10 135 





102 


047 


769 


395 


298 


233 


145 


10 135 


14 


-0 


124 


779 


572 


617 


585 


209 


351 


1Q 137 


-0 


124 


779 


573 


343 


562 


971 


244 


10 137 


15 





154 


446 


315 


576 


951 


606 


777 


10 139 





154 


446 


318 


044 


174 


985 


350 


10 139 



becomes large. 

The problems connected with the computation of ip^k i z ) can largely be avoided if one 
only tries to compute a prediction 



(n) (n,2k) 

9 2k = % 



(43) 



for the first term 7«+2fe+i not used for the construction of either e^} or ip^l\z). For that 
purpose, we have only to set z — in d38|). This yields the following recursive scheme (Eq. 
(3.17) of 0): 



(n) 

9o 

(n) 
92 



(») 
92k+2 



, n e N 
[ln+2] 2 



7«+i 

J«+2) 
92k 



ken. 



n e N 
[7n+2fc+2 — 



9 { 2l +1) ] 2 



7n+2fc+l 

n e N . 



(n) 
92k 



n+2k+2 



-92l +1) ] 2 



7n+2fe+l 



(n+2) 
92k-2 



(44a) 
(44b) 

(44c) 



The main advantage of this recursive scheme over the recursive scheme (J3S|) , from which it 
was derived, is that it only involves numbers but no symbolic expressions. 

In Table |5| we present selected results for the Padc predictions of the first coefficients 
of the perturbation expansion ( p8| ) for AEq{\ 2 ) and ( p9| ) for A£o(/3), which were not used 
in the Pade approximants e 



,(n-2[n/2I,2|In/2]) 



(«-2p/2]) 
2[n/2] 



= -PJ!./S" /21 for 2 < n < 191. The first predictions 



to 6 n +2 and B t 



In/2] 
(n-2[n/2],2[[n/2]) 




to B n +2 were computed with the help of the 
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Table 5: Relative errors R n and lZ n defined in ((45|) and (Q) of 
the Pade predictions for the first coefficients of the perturbation 
series fl28| ) and ( p9| ) not used for the construction of the rational 
approximants. 



n 


Rn 






n n 






2 


-0 295 410 699 






-0 316 117 394 






3 


-0.207 610 910 






-0.218 823 682 






4 


-0.759 683 860 


io- 


1 


-0.833 341 229 


10 


-1 


5 


-0.483 909 816 


10" 


1 


-0.522 231 970 


10 


- 1 


G 


-0.197 254 000 


10" 


1 


-0.218 362 310 


10 


- 1 


7 


-0.120 176 754 


10" 


1 


-0.130 908 858 


10 


- 1 


50 


-0.258 379 657 


10" 


14 


-0.273 374 025 


10' 


-14 


51 


-0.134 007 443 


io- 


14 


-0.141 313 446 


10 


-14 


52 


-0.658 949 507 


io- 


15 


-0.696 514 376 


10 


-15 


53 


-0.341 294 856 


io- 


15 


-0.359 600 798 


10 


-15 


54 


-0.167 932 626 


io- 


15 


-0.177 342 673 


10 


-15 


1 on 


fl ^97 P,7A 71 7 


i n~ 


29 


341 AAQ 


1 n 


-29 


101 


-0.166 900 610 


io- 


29 


-0.173 636 626 


10 


-29 


102 


-0.827 630 688 


io- 


30 


-0.862 097 827 


10 


-30 


103 


-0.421 400 183 


io- 


30 


-0.438 256 801 


10 


-30 


104 


-0.208 999 229 


io- 


30 


-0.217 623 211 


10 


-30 


140 


-0.354 821 178 


10" 


-41 


-0.367 523 224 


10 


-41 


141 


-0.179 775 658 


io- 


-41 


-0.186 000 741 


10 


-41 


142 


-0.893 606 384 


io- 


42 


-0.925 378 776 


10 


-42 


143 


-0.452 675 070 


io- 


42 


-0.468 247 723 


10 


-42 


144 


-0.225 028 932 


io- 


42 


-0.232 976 269 


10 


-42 


187 


-0.295 089 701 


10" 


-55 


-0.304 045 084 


10' 


-55 


188 


-0.146 895 157 


io- 


-55 


-0.151 456 461 


10 


-55 


189 


-0.741 750 230 


io- 


-56 


-0.764 148 505 


10 


-56 


190 


-0.369 259 543 


io- 


56 


-0.380 666 936 


10 


-56 


191 


-0.186 438 979 


io- 


56 


-0.192 040 968 


10 


-56 



recursive scheme (|44|). For the sake of readability, we present in Table ^ only the relative 



errors 



and 



Rn = * 7 =^ (45) 

On+2 



K n = ^ . (46) 

£>n+2 

If the input data of the recursive scheme ( f44| ) are exact rational numbers as the coefficients 
b n or £>„, then the predictions can be computed free of rounding errors. However, it turned 
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out that that the predictions computed in this way were huge rational numbers which slowed 
down computation considerably. Therefore, we used the floating-point arithmetics of Maple 
with an accuracy of 600 decimal digits for the computation of results presented in Table |^. 

The results in Table |f| show that the first coefficients not used for the construction of the 
Pade approximants e^jwij ? ^ = -Fj^/jjj ^ can be predicted with remarkable accuracy by 
the recursive scheme (|44|) if n is sufficiently large. Moreover, the agreement of the relative 
errors R n and lZ n is stunning. However, for our purposes most important is the observation 
that all relative errors in Table | are negative which is in agreement with the Stieltjes 
inequality (|36|). 

Thus, the prediction results of this section show that there is no significant difference 
between the prediction of coefficients b n of the perturbation series (p8|), whose Stieltjes 
nature we want to establish, and the prediction of the coefficients B n of the analogous quar- 
tic anharmonic oscillator perturbation series (|2^), whose Stieltjes nature was established 
rigorously by Simon (Theorem IV. 2.1 of |2|]). Accordingly, these results provide further nu- 
merical evidence that the perturbation series (|2^) for the ground-state energy shift AEo(X 2 ) 
is indeed a Stieltjes series. 



5 Summary and Conclusions 

As shown by countless monographs or articles, the mathematical theory of the conventional 
Hermitian Hamiltonians of quantum mechanics is well established and has reached a high 
degree of sophistication. Moreover, the divergence and the summation of the perturbation 
expansions resulting from these Hamiltonians is also comparatively well understood, in par- 
ticular in the case of the Pade summation of Stieltjes series. In contrast, the rigorous math- 
ematical theory of non- Hermitian, 'PT-symmetric Hamiltonians is virtually nonexistent. It 
is not known whether divergent perturbation expansions resulting from such Hamiltonians 
can be summed to yield a uniquely determined result. Consequently, the best we can do 
for the moment is to perform numerical studies from which we can try to draw general 
qualitative conclusions. 

The main intention of this paper is to provide numerical evidence that the perturbation 
series ( p8| ) for the ground-state energy shift AEq(X 2 ) is a Stieltjes series, because this would 
guarantee that certain subsequences of the Pade table constructed from the partial sums of 
this perturbation series converge to a uniquely determined summation result, as discussed 
in Section ||. 

If the Pade approximants are computed with the help of Wynn's recursive algorithm ( |l3| ) 
according to (|l^) and (|l7|) - as it is done in this paper - and if the series to be transformed 
is the moment expansion of a Stieltjes function, then the Pade approximants must satisfy 
inequality (pTj). As shown in Table |^, the Pade summation results for the perturbation series 
( p8| ) for the ground-state energy shift AE (X 2 ) as well as for the analogous perturbation 
series (|2^) for the ground-state energy shift A£o(/3) of the quartic anharmonic oscillator, 
which is known to be a Stieltjes series, satisfy this inequality. Moreover, the divergence of 
the two perturbation expansions (^8|) and (^9|) as well as as the convergence of their Pade 
summation results is virtually indistinguishable if the two coupling constants A and (3 satisfy 
(i- 

If a Pade approximant Pm{z) constructed from the partial sums (|22|) of the moment 
expansion for a Stieltjes function F(z) is expanded in power series around z = according 
to ([35]), then this series must be strictly alternating for z > just like the Stieltjes series for 
F(z). Moreover, the coefficients \$J m ^ of this expansion are bounded in magnitude by the 
Stieltjes moments /i„ according to the inequality (|36|). 

Thus, via the inequality ( |36"| ) it can be checked whether a moment expansion of the 
type of (|2(]) is a Stieltjes series. However, the accuracy-through-order relationship ( [l2] ) 
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implies that this inequality is by default satisfied for all indices n < I + m, and only for 
n > I + m + 1 do we obtain useful information. In particular for large values of I and m, the 
symbolic construction of Pade approximants and their subsequent expansion may become 
quite demanding both with respect to time and memory. 

These computational problems can to some extent be overcome by expressing the Pade 
approximants P^ +n (z), that can be computed with the help of Wynn's epsilon algorithm, 
according to ( |37j ) by the partial sum from which it was constructed plus the transformation 
term z n+2k+1 ip^} (z). The quantities tp$ (z) can be computed recursively with the help 
of (|3^), and their computation is less demanding than the recursive computation of the 
corresponding Pade approximants P^ +n (z). Moreover, the Taylor expansion of ip^l\z) 

according to (fltf ) yields the desired coefficients ^n+2k+u+i w ^ n v — 0- 

In Tables|3| and ||, we proceed as described above and confirm the validity of the Stieltjes 
inequality d36| ) in the case of the first 15 coefficients of the quantities ygg^(A 2 ) and <£>g6 (/?) 
corresponding to the Pade approximants P33 (A 2 ) and P33 (/3), respectively. In this context, 
it is remarkable that the quality of the predictions is virtually indistinguishable although 
the coefficients b n of the PT-symmetric perturbation series grow much more rapidly than 
the coefficients £>„ for the quartic anharmonic oscillator. 

For larger values of n and k the symbolic computation of the quantities tp^kiz) via 
the recursive scheme ( E}8[ ) becomes quite demanding. In such a case, it is much simpler 
to compute only the prediction for the first series coefficient not used for the computation 
of ip^(z). This can be done with the help of the recursive scheme (|44|). In Table [5] we 
show that all first predictions, which we can compute from the coefficients b v and B v with 
1 < v < 192, satisfy the Stieltjes inequality (|36|). Moreover, the quality of the prediction 
results is again virtually identical. 

We are of course aware that numerical results cannot replace rigorous mathematical 
proofs. Nevertheless, we believe that our numerical experiments are both interesting and 
useful, and that they provide considerable evidence that the perturbation series ( p8| ) for 
the ground-state energy shift AEq(X 2 ) of the PT-symmetric Hamiltonian ([!]) is indeed a 
Stieltjes series, which would imply that the Pade summation of this divergent perturbation 
series converges. 
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